An interface phase transition induced by a driven line in 2D 
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The effect of a localized drive on the steady state of an interface separating two phases in coex- 
istence is studied. This is done using a spin conserving kinetic Ising model on a two dimensional 
lattice with cylindrical boundary conditions, where a drive is applied along a single ring on which 
the interface separating the two phases is centered. The drive is found to induce an interface spon- 
taneous symmetry breaking whereby the magnetization of the driven ring becomes non-zero. The 
width of the interface becomes finite and its fluctuations around the driven ring are non-symmetric. 
The dynamical origin of these properties is analyzed in an adiabatic limit which allows the evaluation 
of the large deviation function of the driven-ring magnetization. 
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The effect of local drive on the properties of an inter- 
face separating two coexisting phases has recently been 
explored as a simple example of systems driven away from 
equilibrium. Much of the attention is due to the surpris- 
ing experimental results on colloidal gas-liquid interface 
subjected to a shear flow parallel to the interface [1 . It 
was found that the shear drive applied away from the 
interface, strongly suppresses the fluctuations of the in- 
terface, making it smoother. This long-distance effect 
of the drive is due to long-range correlations that char- 
acterize driven systems [2]-[8]. An interesting theoretical 
approach for studying this phenomenon has been intro- 
duced by Smith et al. who considered a two dimensional 
version of the system, and modeled it by an Ising lattice- 
gas below its transition temperature [9] . Using spin con- 
serving Kawasaki dynamics and applying shear flow at 
the boundaries parallel to the interface, it was observed 
that the interface indeed becomes narrower although its 
width still increases with the length of the interface. In 
closely related works, the effect induced by a current car- 
rying line on a neighboring non-driven one has also been 
analyzed [TqHT3]. 

In this Letter we consider a drive localized along an in- 
terface which separates two coexisting phases, and study 
the resulting interface properties. This is done using a 
two dimensional Ising model on a square lattice with 
cylindrical boundary conditions (Fig. [T]), that evolves 
under spin conserving dynamics. The drive acts along 
the ring around which the interface is centered. We 
find that the drive induces an interface phase transition 
which involves spontaneous symmetry breaking, resulting 
in a non-zero magnetization of the driven ring. In this 
transition, the macroscopic 2D steady state remains un- 
changed, however spontaneous symmetry breaking takes 
place involving the steady state of a ID stripe centered 
on the driven ring. This is in sharp contrast with an 
equilibrium setup of an interface subjected to a localiz- 
ing potential along a ring, where the ring magnetization 
vanishes at all temperatures, and no interface sponta- 
neous symmetry breaking takes place. Moreover, we find 
that the drive suppresses the fiuctuations of the inter- 
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FIG. 1. (color online) Square lattice with cylindrical bound- 
ary condition, with the drive on the central ring and the 
boundary conditions indicated. 



face, leading to an interface with a finite width which 
does not scale with the system size. Also, due to the 
broken symmetry on the driven ring, the interface fluc- 
tuations are highly asymmetric. The interface fluctuates 
more strongly into the bulk phase whose magnetization 
is opposite to that of the driven ring. These results are 
first demonstrated by numerical simulations. The model 
is then analyzed in a special limit which allows analytical 
computation of the large deviation function (LDF) J4] 
of the magnetization of the driven ring, demonstrating 
the existence of the spontaneous symmetry breaking. 

To proceed, we consider Ising spins cr = {a^} on sites 
r = (x, y) of an L X (2M + 1) square lattice, with periodic 
boundary condition in the x-direction while the two open 
boundaries, y = ±M, are coupled to rows from above 
{y = M and below (?/ = — M — 1), respectively, with 
fixed spins: crx^±M±i = Tl (Fig. [T]). The model has 
nearest neighbor ferromagnetic interactions and a drive 
is introduced by a force field E = {E^ 0) applied on the 
y = ring. The field favors the positive spins to move 
counter clockwise along the ring, and as a result drives 
the system out of equilibrium. 

There is more than one way to incorporate the drive in 
the dynamics, and unlike a dynamics satisfying detailed 
balance, the steady state depends on the precise choice 
of the rates We choose a modified Metropolis algo- 
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around the driven ring (denoted by a dashed hne) for L = / 
2M = 100 at T = 0.85Tc. Blue (dark) and yehow (hght) 

squares indicate + and - spins, respectively. 3 ^^^.^^^ ^.^^ evolution of mo for L = 2M = 100 at 

T — 0.6Tc. Time is measured in Monte Carlo steps. 



rithm [3] where in every step a pair of nearest neighbor 
sites r and are chosen at random and their spins are 
exchanged with probability min{l, exp(— /^A'H)}, where 
P is an inverse temperature and AH is the energy differ- 
ence between the final and initial configurations. Thus, 
for exchanging cFy and CTr', 

^ I Ai7 - (dr - (Jr')(v' - r) • E if r,r' G 0th ring, 
I Ai^ elsewhere, 

where Ai^ is calculated using the Ising Hamiltonian H = 
~J^{yy')^y^y'^ with J > 0. One Monte Carlo time 
step is constituted of L(4M + 1) such updates. In all the 
numerical results presented in this Letter we use a large 
driving field E > lOJ. 

In absence of a driving field, the model is in equi- 
librium. At sub-critical temperatures {T < Tc — 
2.2692J//cb), the equilibrium state is composed of two 
oppositely magnetized phases, separated by an inter- 
face. For an initial configuration with zero overall mag- 
netization, the magnetization profile in the y direction, 
rriy = ^/L^^dx^y^ is antisymmetric with respect to 
y = 0. The interface fluctuates symmetrically around 
the driven ring, leading to zero magnetization on the ring, 
mo = 0. In the large L, M limit with fixed aspect ratio 
L/M^ the width of the interface scales as a/Z [T6 . 

Introducing a drive does not modify the overall macro- 
scopic structure of the steady state. As is naively ex- 
pected, the steady state is still composed of two oppo- 
sitely magnetized phases separated by a fluctuating in- 
terface around y = 0. However, numerical studies of 
the model reveal some profound changes in the struc- 
ture of the interface itself. In particular we find that 
(a) in the thermodynamic limit the magnetization of the 
driven line, mo, is non-zero, taking one of two oppositely 
directed values. It thus breaks the <Jr — <j_r symmetry 
of the model, (b) The interface is localized around the 
driven line and its width stays finite in the thermody- 
namic limit, and (c) the fluctuations of the interface into 
the two bulk phases are highly asymmetric, with more 
pronounced fluctuations into the phase whose magneti- 
zation is oppositely directed to that of the driven line. 
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FIG. 4. Average time between consecutive switches ts{L) for 
2M = L, and T = 0.6Tc 



In Fig. |2]we present two typical microscopic configu- 
rations of the model. It is clearly seen that the driven 
line is predominantly occupied by either positive or neg- 
ative spins representing its two possible ordered states. 
As the system evolves, the magnetization mo fluctuates 
around one of the non-zero values for a long time. It then 
switches to the oppositely magnetized state over a much 
shorter time scale, as shown in Fig. [Sj 

The numerical result for the average time between two 
successive such switches, tg, are shown in Fig. [4] The 
data suggest that tg grows exponentially with L, with 
ts ~ exp(0.06L). The data for each L is averaged over n 
number of switches that are observed in available com- 
putation time (n varies from around 12000 to 10 as L 
changes from 30 to 130); n decreases with L, yielding in- 
creasing error bars of order with L. Although the 
range of the system size studied is insufficient for a con- 
clusive evidence of an exponential growth, this form is 
justified by the theoretical results presented below. The 
exponential growth implies that in the thermodynamic 
limit, the two non-zero values of mo correspond to two 
thermodynamically stable phases. 

The width of the interface is evaluated by averaging 
\y\ weighted by the derivative drriy/dy that peaks at the 
interface position. The result is shown in the inset of Fig. 
|5]for both driven and non-driven case. A comparison of 
the two cases clearly indicates that the interface fluctu- 
ations are drastically reduced in the presence of drive. 
As will be shown by the theoretical analysis presented 
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FIG. 5. (color online) The average magnetization profile my 
corresponding to the two phases, close to ^ = 0. The asymme- 
try around y = is clearly seen. The profiles are generated on 
a 100 X 101 lattice at T = 0.85Tc, averaging over 10^ config- 
urations at regular intervals of 1000 Monte Carlo steps. The 
figure in the inset shows the growth of the width of the inter- 
face with increasing length L for zero drive (top curve) and 
for driving strength E = 10 J (lower curve). 



below, the width of the interface remains finite at large 
L. Such smoothening of the interface has also been ob- 
served in presence of global drive parallel to the interface 
[17] . The interesting difference here is that the interfacial 
fluctuations are asymmetric, resulting in an asymmetric 
magnetization profile around the driven line (see Fig. [5| . 

In order to make an analytical analysis of the model 
feasible, we generalize the model by introducing a param- 
eter 7 that controls the dynamical rate of the processes 
involving spin exchange between the y = ring and the 
neighboring rings y = ±1. For these processes the rate 
becomes min{7, 7exp (— /3AH)}, with 7 > 0. The other 
rates remain unchanged. This does not modify the steady 
state of the equilibrium case {E = 0) but it helps analyz- 
ing the non-equilibrium steady state. We now consider 
the steady state in the following special limit: (a) slow 
exchange rates (7 <C L~^) between the driven and the 
neighboring rings, (b) an infinite driving field {E 00), 
and (c) low temperature (exp(— /3J) <C 1). We show 
below that in this limit the stationary probability distri- 
bution P (mo) of the magnetization mo of the driven line 
has the form P(mo) cx exp(— L(/)(mo)). The large devi- 
ation function, 0(mo), is then computed and shown to 
possess two degenerate minima at non-vanishing values 
of the magnetization mo = imj (see Fig. [6|, implying 
a spontaneous symmetry breaking on the ring. In ad- 
dition, the LDF yields an exponential flipping time in 
L between positive and negative magnetization for finite 
systems due to the finite barrier between the two minima. 

We proceed by noting that due to the slow exchange 
rate 7, there are no significant exchanges between the 
driven line and its neighboring rings on a time scale 
t ^ At = {jL)~^. On such time scale the lattice may be 
considered as composed of three subsystems: the driven 
line, and the upper u {y > 0) and lower i {y < 0) sub- 
lattices. They evolve while keeping their own specific 



FIG. 6. The large deviation function (j){mo) calculated using 
Eq. m and ([6|. 

magnetization mo, m^, and m^ unchanged, reaching the 
steady state corresponding to fixed subsystem magnetiza- 
tion. On a longer timescale, t ^ At, the magnetizations 
mo, m^i, and mi evolve as spins are exchanged between 
the subsystems. 

We now define a coarse-grained time variable r = t/ At 
such that the subsystem magnetization evolves with in- 
creasing r, however at any given r each subsystem is 
effectively in the steady state corresponding to its mag- 
netization. This separation of slow and fast processes 
is analogous to the adiabatic approximation in quantum 
mechanics [18^ , and has also been applied in related mod- 
els [19l[20]. 

Let us characterize the steady states corresponding to 
fixed sub-system magnetization mo, rriu and m^. First 
consider the driven ring. In the limit E ^ 00 the dy- 
namics within this ring is independent of the two other 
subsystems, and reduces to that of the Totally Asym- 
metric Simple Exclusion Process (TASEP). In its steady 
state all spin configurations with fixed magnetization mo 
are equally probable, leading to uniform magnetization 
and zero spin- spin correlation along the driven line. This 
steady state is reached in a time of 0(L^/^) which, for 
7 <C L~^^ is smaller than the typical time of exchange 
processes between the driven ring and its neighboring 
ones. Then, the driven ring provides an effective bound- 
ary magnetic field Jmo on the u and £ subsystems. Thus 
the steady state of these two subsystems is the equilib- 
rium state of the Ising model subjected to a boundary 
field. The boundary field results a magnetization profile, 
m^, which for large \y\ approaches the bulk magnetiza- 
tion values —rriB and ttib for the u and £ subsystems, 
respectively. The length scale of this approach is of the 
order of the spin-spin correlation length ^(T) of the 2D 
Ising model. Since this length is finite at all tempera- 
tures except at Tc this demonstrates that the width of 
the interface remains finite for large L. 

Let Pr{mo)dr be the probability of the driven line 
magnetization to have value mo between coarse-grained 
time r and r + (ir, while rriu and m^ have already reached 
stationary values. The probability function evolves as 
spins are exchanged between the subsystems. At each 
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exchange process between the driven ring and the bulk, 
mo changes by ±21 L. Let p{mo) and q{mo) be the in- 
creasing and decreasing rates of mo, respectively. Then 
the dynamics of mo is that of a random walker with posi- 
tion dependent forward and backward jump rates p{mo) 
and q{mo), respectively, and with boundary condition 
P^(mo) for |mo| > 1. 

The stationary distribution of this motion is an equi- 
librium distribution function. 



P(mo) = P(0)exp[-L(/)(mo)]. 



(1) 



The LDF ^(mo) is an even function of mo and it 
can be determined using the detailed balance condition 
p{mo)P{mo) = q{mo + 2/L)P{mo + 2/L), which for 
< mo < 1 yields 



P mo 



2n\ 



1)] 



k=l 



(2) 



with n = l,...,L/2. In the large L limit this yields the 
LDF for m > 0, 



^(mo) 



1 



dm In 



p (m) 



q{m) 



(3) 



with (/)(— mo) = (/)(mo). 

The rates p{mo) and q{mo) are determined as fol- 
lows: consider a spin exchange process between the 
driven ring and its two neighboring ones, in which 
the microscopic configuration changes from a to a' 
and mo increases by 2/L. The rate of this process 
is 00 {a ^ cr^) P (tT|mo, — m^, ttib) where uo {cr ^ cr') = 
min{l, exp(— /3AH)} is the Metropolis success rate in 
coarse grain time variable r and P {a\mo, —mB^rriB) is 
the steady state probability of configuration a corre- 
sponding to subsystem magnetization mo, m^ = — m^ 
and m£ = m^. Summing over all such exchanges one 
obtains 

p (mo) = ^uj{a ^ a')P {cr\mo, -rriB.mB) , (4) 



where the sum is over configurations cr' whose mo is 
higher than that of a by 2/L. The magnetization de- 
creasing rate g(mo) is readily obtained by noting that 
due to the invariance of the dynamics to space-time in- 
version, (Jr — cT-r, one has q{mo) =p{—mo). 

In the slow exchange limit 7 <C L~^^ the probability 
P (o-|mo, — m^, tub) can be expressed in terms of proba- 
bility of the subsystem configurations as 

P(o-|mo, -TUB, tub) c^P{cro\mo)P (t7^| - m^, mo) 

P {ai\mB,mo) , (5) 

where cro, cr^ and cr^ are the microscopic spin config- 
urations of the three subsystems corresponding to the 
configuration a. Here, P((To|mo) is the steady state 



+ 1 + 1 + 
+ + + 



+ |-| + 
+ + + 



FIG. 7. (color online) Spin exchange process which yields 
the leading order contribution to p(mo). The driven ring is 
denoted by a dotted line, and the pair of spins exchanged are 
indicated. 



distribution of the driven line with fixed magnetization 
mo, which is the same as the steady state of a TASEP, 
and P {(Tul —rriB^mo) and P(cr^|mB,mo) are the equi- 
librium distribution of the other two subsystems. 

In general, calculating all the terms in Eq. Q is not 
straightforward. However, the calculation becomes feasi- 
ble in the low T limit where these rates may be expanded 
in powers of exp (— /3J). 

In order to keep track of the terms in this expansion 
it is convenient to generalize the model by considering 
an interaction strength between the driven ring and 
its neighboring ones as Ji < J. It is easy to see that 
the leading contribution to p{mo) in Eq. Q results 
from the exchange process shown in Fig. [7| where both 
subsystems u and £ are in their respective ground state, 
m£ = —rriu = 1. For this process uj {a ^ a') = 
exp(— 2/3(J + Ji)) and P {a\mo^ —mB^rriB) = 
(1 + mo)^ (1 - mo) /sj [1-0 (e-^^^)] . Higher or- 
der contributions can be determined similarly from other 
exchange events. Computing p(mo) up to (9(exp(— 6/3J)) 
yields, for —1 < mo < 1, 



P (mo) = + ^0)' (1 - mo) e-'^^^ 

O - 

1 r 



-2/3J 



(1 + mo)' (1 - mo) (2e-2^^i^° + ^^^Jimo^ 
+ 2 (1 + mo) (1 - mo)' {e'^^^' + e'^^^^^) 

+ (1 - mo)' 6^^^^°] e-^^^ + O {e-'^') . (6) 

The LDF (j){mo) calculated using the rate in Eq. ^ is 
plotted in Fig. [6] for /3Ji = /3J = 3/4 (T O.GTJ. This 
function has two minima which correspond to the two 
thermodynamic phases with non-zero mo- The average 
time Ts for the magnetization to jump from one minimum 
to the other is proportional to the exponential of the bar- 
rier height between them. In terms of Monte Carlo steps 
this switching time tg ~ 7~^Ts ~ 7~^exp(eL) where e 
is the barrier height. For the parameters of Fig. [6] one 
has e = 0.18 which is of the same order as that obtained 
numerically in Fig. [4j For a better comparison higher or- 
der terms in the low temperature expansion are required. 
The asymmetry in the fluctuations of the interface and 
the magnetization profile in Fig. [5] is a consequence of 
the non-zero values of mo- 
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The analysis presented in this Letter demonstrates 
that a local drive can induce a phase transition which 
involves spontaneous symmetry breaking of an interface 
separating two coexisting phases. It would be interesting 
to consider other boundary conditions which would allow 
the interface to detach from the driven ring. This would 
correspond, for example, to studying the model with pe- 
riodic boundary conditions in both the x and y directions. 
In this case the model exhibits two interfaces, and pre- 
liminary studies have shown that either one of them is 
attracted by the driven ring resulting in a macroscopic 
symmetry breaking, in addition to that of the interface 
[21]. This will be addressed in a future publication. 
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